How do I compute the Lyapunov exponents for Lorenz system

I want to find all the max or min values in a range, from a numerical function [t,x(t)], numerical solution of a differential equation.

For example if the solution was periodic in time, I want to find the max values of x in a certain intervalot t